Controlling polarization at insulating surfaces: 
quasiparticle calculations for molecules adsorbed on insulator films 
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By means of quasiparticle-energy calculations in the Go Wo approach, we show for the prototypical 
insulator/semiconductor system NaCl/Ge(001) that polarization effects at the interfaces noticeably 
affect the excitation spectrum of molecules adsorbed on the surface of the NaCl films. The magnitude 
of the effect can be controlled by varying the thickness of the film, offering new opportunities 
for tuning electronic excitations in e.g. molecular electronics or quantum transport. Polarization 
effects are visible even for the excitation spectrum of the NaCl films themselves, which has important 
implications for the interpretation of surface science experiments for the characterization of insulator 
surfaces. 
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On the nanoscalc, materials often reveal extraordinary 
features. To harness this potential it is essential to grow 
or manufacture nanostructures in a controlled way. Ul- 
trathin insulating films are an example for which this has 
been achieved on metal and semiconductor surfaces. We 
have recently shown that these films develop new and 
unique properties as their thinness approaches the limit 
of a few atomic layers and that such supported ultrathin 
films should be regarded as new nanosystems in their own 
rights Here we go one step further and demonstrate 
by means of first principles calculations that control over 
the film thickness means control over the polarization of 
the film. This in turn gives access to properties on the 
film's surface, for example the energy levels of molecu- 
lar adsorbates, which are relevant in the context of e.g. 
catalysis, molecular electronics, or quantum transport. 

The fact that ultrathin insulator films offer a new per- 
spective of control at the nanoscale is increasingly be- 
ing recognized. Repp et al, for example, have recently 
demonstrated that gold atoms adsorbed on NaCl/Cu can 
be reversibly switched between the neutral and negative 
charge state Q. Alternatively, the charge state of Au 
atoms on MgO films can be controlled by the film thick- 
ness 3 . For planar molecules adsorbed on Cu-supported 
NaCl films, the molecular orbitals can be resolved spa- 
tially and energetically [H, 0| by scanning tunneling mi- 
croscopy (STM) and spectroscopy (STS), and even re- 
actions can be followed Interestingly, STS exper- 
iments performed on pentacene molecules adsorbed on 
NaCl/Cu(lll) show a significant influence of the film 
thickness on the molecular gap [f|. 

Ultrathin insulator films have also developed into 
highly valuable and intensively studied model systems 
for characterizing insulating surfaces. The study of in- 
sulator surfaces has proven difficult due to their lack 
of conductivity, which severely limits the range of ap- 
plicable surface science techniques. Ultrathin insulator 
films grown on conducting substrates offer a solution to 



FIG. 1: (Color online). Structure of NaCl films (2-4 layers) 
on Ge(00f) with a CO adsorbate. The spheres denote the 
atoms: Ge (medium magenta/gray), Na (large blue/dark), 
CI (large green/gray), C (small black), O (small red/gray). 



this dilemma because the films can exchange electrons 
with the substrate by tunneling 0, @|. Caution has to 
be applied, however, when transferring thin-film results 
to the surfaces of technological interest. The properties 
of ultrathin films may deviate considerably from those of 
macroscopic films |l| and the excitation spectrum may 
be affected by the polarization effects presented here. 

In this Letter we will address both of these points by 
means of GqWq quasiparticle-energy calculations [l(| 
for the example of a prototypical insulator-semiconductor 
interface (NaCl/Gc(001)) and CO as a model adsor- 
bate. Supported NaCl films are well-behaved model sys- 
tems for studying the properties of insulator surfaces 
@, El El El- Although they arc mostly grown on 
metals, notably Cu 0.G3, Ge(OOl) is the substrate of 
choice for stud ying insulator/semiconductor interfaces 
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12 . 1 14 , Il5l . Ilq . In recent years, these films have 
also attracted increasing interested in the context of STM 
and STS studies of atomic and molecular adsorbates 

We use density-functional theory (DFT) in the local- 
density approximation (LDA) to determine the atomic 
structure of NaCl films on Ge(OOl). The electronic exci- 
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tation spectrum is calculated with many-body perturba- 
tion theory in the GW approach as perturbation to the 
LDA ground state (henceforth denoted G W @LDA). 
The GW approach has not only become the method of 
choice for calculating quasiparticle excitations in solids 
as probed in STS or direct and inverse photoe- 
mission, but also includes long-range polarization effects. 
In a bulk material these encompass the screening of ad- 
ditional charge. At a surface or an interface, however, 
the abrupt change in dielectric constant gives rise to a 
net build up of charge, so called image charges. This net 
polarization acts back on the additional charge and in- 
creases in strength with decreasing distance between the 
additional charge and the interface. These polarization 
effects are absent from the most common density func- 
tionals (such as the local-density or generalized gradi- 
ent approximation, exact-exchange and hybrid function- 
al), but enter the GW self-energy (E=iGW) through 
the screened Coulomb potential W . The application of 
the GW method is therefore necessary to capture these 
effects [II El El, [II . 



Polarization or image effects are present at any sur- 
face or interface, but are most commonly associated 
with metal surfaces, where the ratio in dielectric con- 
stants is largest. In supported ultrathin films, however, 
a charged excitation (e.g. an electron added to the 2n* 
level of CO on NaCl/Gc(001)) polarizes two interfaces 
(vacuum/NaCl and NaCl/Ge(001)). The combination of 
dielectric constants and film thickness therefore controls 
the strength of the polarization effects from the semicon- 
ductor/insulator interface to the insulator surface. 

Before we address the polarization effects at the 
NaCl/Gc interface in detail we briefly describe its atomic 
structure, which had not been determined previously. 
The DFT-LDA calculations were performed with the 
SFHIngX code [2lj]. We employ a plane- wave basis set 
(40 Ry cutoff) and norm-conserving pseudopotcntials. 
The NaCl/Ge system is modeled in the repeated slab- 
approach with a 6-layer Ge slab at the experimental lat- 
tice constant (saturated by hydrogen atoms on the bot- 
tom side) plus a varying number NaCl layers on the top 
side (22)1. Increasing the Ge thickness to 12 layers pro- 
duces no significant change in the atomic or electronic 
structure of the NaCl films. In agreement with experi- 
mental indications [l6| , we find that the Ge dimers of the 
clean Ge(001) surface prevail below the NaCl film, giving 
rise to a 2 x 1 surface lattice [H| . The Ge dimers below 
the NaCl film remain asymmetric, but with a smaller 
tilting angle (10°) compared to the free surface (19°). 
The adhesion of the film is dominated by the electro- 
static interaction between the ions in the film and the 
partial charges developing at the buckled-dimer surface 
of Ge(001). The relaxation pattern in the bottom layer 
follows the electrostatic profile of the Ge surface, which 
attracts the ions next to the dimer, but repels those above 
the inter-dimer troughs (cf. Fig. [T]). The corrugation 
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FIG. 2: a) 5<t-2tt* splitting of CO on NaCl/Ge(001) at the 
LDA and GoW-oQLDA levels of theory, and b) GqWq correc- 
tions to the 5a and 2-7T* molecular orbitals as a function of 
the NaCl film thickness. 



in the higher layers is induced by the bottom layer and 
quickly flattens out as the thickness increases. 

Next, we use the adsorption of a small molecule to 
probe the effect of the interface polarization at the film's 
surface. For this purpose, we placed a single CO molecule 
in the 2x1 surface unit cell (we have no indications 
for relevant changes at lower coverages). CO physisorbs 
perpendicular to the NaCl(OOl) surface with the C-end 
down (cf. Fig. [J). The CO axis tilts along (110), i.e. 
the Ge dimer, by 2.6° for 2 monolayer(ML) NaCl, 0.7° 
(3 ML), and 0.1° (4 ML), respectively. The adsorption 
energy of 0.28 eV is the same for the two incquivalent Na 
sites at the surface and does not depend on the film's 
thickness to within 0.01 eV. It also agrees with the value 
for a thick free-standing NaCl slab. The molecular states 
give rise to flat bands in the band structure and do not 
hybridize with NaCl states. In the following, we focus on 
the molecular gap (given by the ha-2n* splitting) at the 
r point. 

The Go Wo calculations for the electron addition and 
removal spectra were performed with the gwst code 
[23L [241 . [23 ■ For the correlation (exchange) self-energy, 
a 14 Ry (28 Ry) plane-wave cutoff and a 6 x 3 x 1 It- 
point sampling was used. State summations included 
2500 bands (81 eV above the Fermi level). To correct 
for artificial polarization effects in the repeated-slab ap- 
proach, a "finite vacuum" correction was applied [26|. 

In Fig. [2j we compare the molecular gap of CO com- 
puted for LDA and G W @LDA for films of 2-4 ML 
thickness. At the level of LDA, we observe a small re- 
duction of the molecular gap with increasing film thick- 
ness due to the thickness-dependent structural changes 
in the surface of the NaCl film. The Go Wo corrections, 
however, reverse this trend and introduce a significant 
increase in the quasiparticle gap from 12.48 eV (for 2 
ML) to 12.80 eV (4 ML). The other CO orbitals exhibit 
analagous thickness-dependent shifts (not shown) . These 
gaps are significantly smaller than for CO on a pure 
NaCl surface (13.1 eV for 6 ML NaCl, no Ge) or the free 
molecule (15.1 eV). Similar reductions have been found 
in GqWo calculations for molecules adsorbed on insula- 



3 



tors |27j or semimetals [28| and are a result of surface 
polarization or charge transfer to the molecule pij . 

We will now demonstrate that the reduction of the 
5cr-27r* splitting is caused by polarization effects at the 
two interfaces. Inspection of Fig. [3b) reveals that the 
thickness-dependent changes are of similar magnitude, 
but opposite sign. Such a behavior is characteristic for 
long-range polarization effects [26|, [28j|. To illustrate this, 
we split the electron addition/removal into two steps. 
First, the charged state is created on a free molecule. 
Secondly, we consider the electronic polarization of the 
NaCl/Gc substrate. If the lifetime of the charged state on 
the molecule is long compared to the electronic relaxation 
time of the substrate, we may treat the rearrangement 
of charge with classical electrostatics. The polarization 
lowers the energy of the charged state E N ±i by ^qAV, 
where AV is the polarization-induced change in the elec- 
trostatic potential and q the charge. The energy of the 
hole state (En — En-i) thereby increases whereas that 
of the electron state (En+i — En) becomes smaller. 

To substantiate this picture, we estimate AV outside 
the supported NaCl films. For this, the Ge substrate, the 
NaCl film, and the vacuum region arc replaced by homo- 
geneous dielectric media (e—14, 2.8, and 1, respectively) 
with abrupt interfaces. The thickness of the NaCl region 
is taken to be 2.8 A/layer. The effective polarization is 
then computed using the image-charge method. For AV, 
we then take the value of the image potential at the po- 
sition of the CO molecule outside the surface. Since this 
position is somewhat ambiguous due to the spatial ex- 
tent of the CO orbitals we determine it by requiring that 
the reduction of the 5<r-2ir* splitting of 2 eV at the bare 
NaCl surface is reproduced by the model. This yields a 
value of 1.5 A, and the model then gives 5(t-2tt* splittings 
of 12.6, 12.7, and 12.8 eV for the 2, 3, and 4 ML films, 
respectively, in good agreement with the values from our 
G W @LDA calculations. 

On the experimental side, a gap reduction as a function 
of film thickness has been reported for STS experiments 
on pentacene molecules adsorbed on NaCl/Cu(lll) [B|. 
In STS, the molecular states give rise to tunneling reso- 
nances and the observed tunneling gap amounts to 3.3, 
4.1, and 4.4 eV for NaCl films of 1,2, and 3 ML in thick- 
ness (the gap of pentacene in the gas phase is 5.3 eV). 
The overall trend as well as the magnitude agree very 
well with our G W @LDA calculations for CO/NaCl/Gc 
considering that the ratio of dielectric constants is much 
larger in the NaCl/Cu case. A tunneling resonance gap 
will be hard to observe experimentally for CO/NaCl/Ge, 
however, because the highest occupied molecular orbital 
(5a) is located ~8.5eV below the Ge valence band max- 
imum (vbm) in GoWo@LDA and thus even below the 
NaCl valence band (4-7 eV below the Ge vbm, cf. Fig. [3]). 

Our results highlight that surface and interface po- 
larization effects are important for adsorbed atoms and 
molecules. Supported ultrathin films thereby offer un- 
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FIG. 3: DOS inside the supported NaCl films of 2-4 ML 
thickness at the LDA and the GoWo@LDA level. The shaded 
region indicates the energy position of the Ge bands. The 
valence band maximum of the Ge slab defines the energy zero. 



prcccdentcd opportunities for controlling these effects by 
tailoring both the film's thickness and dielectric constant 
to the desired properties. These additional parameters 
may expedite the design of devices in molecular electron- 
ics or quantum transport, where the distance of a molec- 
ular state to the Fermi level is an important quantity. 

We will now demonstrate that the same polarization ef- 
fects also affect excitations inside the ultrathin film and 
discuss the implications of this for the characterization 
of insulator surfaces. Fig. [3] shows the density of states 
(DOS) projected onto the NaCl film for films with 2-4 
ML for LDA with and without the GqWq corrections. 
The most important change when comparing the LDA 
and G Wo DOS is the shift of the NaCl bands [H rel- 
ative to the Ge states. This is not too surprising, as 
the Go Wo corrections to the bulk band gaps are much 
larger for NaCl (3.3 eV) than for Ge (0.7 eV). Including 
the Go Wo corrections, the top of the film's valence band 
lies ~4.2eV below that of the Ge substrate in excellent 
agreement with ultraviolet photoelectron spectroscopy 
3l| . More remarkable, however, is the change in the 
shape of the NaCl-derived features when going from LDA 
to GoWo@LDA and from 2 to 4 layers. This indicates 
that the Go Wo shifts for the NaCl states are not uniform 
and that corrections derived from bulk calculations are 
not easily transferable to thin films. Excited states are 
instead subject to additional thickness-dependent and 
substrate-specific changes that are neither encompassed 
by a ground-state perspective nor easily derivable from 
bulk Go Wo calculations for the separate fragments alone. 
These thickness-dependent variations should be observ- 
able in high-resolution spectroscopic experiments. 

In LDA the DOS of bulk NaCl (not shown) is already 
attained at a thickness of only three layers. A similar be- 
havior has previously been observed for ultrathin silica, 
hafnia and alumina films [l|. However, this is no longer 
the case when charged excitations (e.g. photocmission or 
tunneling) are treated appropriately. The GqWo@LDA 
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FIG. 4: (Color online) Comparison of the extrapolated Go Wo 
orbital shifts for the 4-layer NaCl film to the image potential 
of the dielectric slab model (cf. text). 



calculations demonstrate clearly that the DOS of ultra- 
thin films differs from that of bulk NaCl (which is iden- 
tical in shape to the LDA DOS for a 4 layer film). This 
implies that caution has to be applied when interpreting 
spectroscopic results. Spectra of ultrathin films are not 
representative of bulk samples. 

The non-uniform GqWq shifts are also a result of the 
interface polarization effects. The position-dependence of 
the self-energy (S = iGW) can be understood with the 
same dielectric slab model introduced above. Solving the 
image-charge model yields the induced potential AV(z) 
as a function of the position in the film z. From this, 
we build the following model self-energy for the occupied 
NaCl states 
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E(z) = AS - -AV(z) 



(1) 



The constant AS encompasses all GW self-energy effects 
that are not due to image effects. Its actual value is not 
important for analyzing the position dependence of the 
self-energy. Based on our GqWq calculations, we estimate 
it to be -1.05 cV for CI 3s and -1.3 cV for CI 3p. In or- 
der to assess the local variation of the self-energy in the 
full GqWq calculation, the valence bands are projected 
onto atomic orbitals. For states that are predominantly 
composed of orbitals from a single atomic layer (>80%) 
we observe a linear dependence between the Go Wo cor- 
rection and the projection weight. By extrapolating 
to 100%, the orbital-dependent quasiparticle corrections 
shown in Fig. [4] are obtained. The variation of the quasi- 
particle shifts throughout the NaCl film (Go Wo DOS in 
Fig. [3]) is reproduced well by a projection-weighted sum 
of these local-orbital contributions. In Fig. [5J we com- 
pare the extrapolated shifts to our model self-energy. The 
good agreement illustrates that the interface polarization 
effects in the Ge/NaCl system are indeed the cause for 
the position dependence in the Go Wo self-energy. 

We acknowledge fruitful discussions with Philipp Eg- 
gert and the Nanoquanta Network of Excellence (NMP4- 
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